Spatial antibunching of photons with parametric down-conversion 
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The theoretical framework behind a recent experiment by Nogueira et al. [Phys. Rev. Lett. 86, 
4009 (2001)] of spatial antibunching in a two-photon state generated by collinear type II parametric 
down-conversion and a birefringent double-slit is presented. The fourth-order quantum correlation 
function is evaluated and shown to violate the classical Schwarz-type inequality, ensuring that the 
field does not have a classical analog. We expect these results to be useful in the rapidly growing 
fields of quantum imaging and quantum information. 
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I. INTRODUCTION 

As current technology advances, more and more atten- 
tion is placed upon Quantum Mechanics to solve future 
problems. Furthermore, quantum systems are capable of 
performing some tasks more efficiently than classical sys- 
tems [0, drawing even more emphasis to quantum tech- 
nologies. In particular, the fields of optical communica- 
tion, optical imaging and optical information processing 
have been appended by the rapidly developing fields of 
quantum communication j|, ||, Q| , quantum imaging || |(| 
and quantum information processing |Q. Thus, the study 
of quantum phenomena promises to be a fruitful enter- 
prise. 

For many years researchers have studied the non- 
classical behavior of light, such as squeezing ||, |jj and 
antibunching |ll|, However, most theoretical and 
experimental investigations deal with time variables only. 
That is, most treatments consider only one spatial mode. 
In a recent review article, Kolobov [fbf demonstrates that 
many quantum phenomena also occur when considering 
spatial variables of the electromagnetic field. Many areas 
of technology stand to benefit from the possible applica- 
tions provided by such quantum phenomena. 

An invaluable tool in these areas of research is the 
generation of entangled photons using parametric down- 
conversion fl4|j . The two-photon state of light exhibits 
non-separable behavior and has been used in 

nearly all quantum information schemes |]l7f . 

Spatial antibunching was recently observed experimen- 
tally by Nogueira et al. [jl8| using spontaneous paramet- 
ric down-conversion (SPDC). In this article, we provide 
a theoretical background for the experiment reported in 
ijbSf . Section || is dedicated to the general introduction of 
temporal and spatial antibunching. In section III we dis- 



cuss the theoretical observation of spatial antibunching 
of photons using a two-photon entangled state produced 
by SPDC, as in p8[. We close with some concluding 



remarks in section IV 



II. PHOTON BUNCHING AND 
ANTIBUNCHING 

It is well known that any state of the electromagnetic 
field that has a classical analog can be described by means 
of a positive nonsingular Glauber- Sudarshan P distribu- 
tion, which has the properties of a classical probability 
functional over an ensemble of coherent states. Because 
of this fact, the normally-ordered intensity correlation 
function for stationary fields must obey the following in- 
equality [pi: 



(T:I(r,t)I(r,t + r) :) < (: I 2 (r, t) :}, 



(1) 



where T : : stands for time and normal ordering. Photon 
density operators are defined as 



I(r,t) = V\r,t)-V(r,i), 



where 



V(r, t) = — 1= ^ Ok^ek^e 

k,o- 



i(k-r-w k t) 



(2) 



(3) 



Ok,cr is the annihilation operator for the mode with wave 
vector k and polarization <r, e^^ is the unit polarization 
vector, n is the quantization volume and u> — ck. 

Expression (Q) is commonly written in the shorter form 



where 



G^iv^v^KG^iy^O), (4) 



G^ 2 \v l ,v 2 ,T) = (T:i(v,t)i{v,t + T):). (5) 
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Since the delayed photon coincidence detection proba- 
bility ■P(ri,r 2 ,r) is proportional to G( 2,2 )(ri, r 2 , r) p)] , 
inequality (^) means that for the class of fields consid- 
ered above, photons are detected either bunched or ran- 
domly distributed in time. Photon antibunching in time, 
characterized by the violation of ([!]), was predicted by 
Carmichael and Walls [[[o| , Kimble and Mandel [[ll] , and 
was first observed by Kimble, Dagenais and Mandel in 
resonance fluorescence [O] . 
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In the space domain, the concept analogous to sta- 
tionarity is homogeneity. For a homogeneous field, the 
expectation value of any quantity that is a function of po- 
sition is invariant under translation of the origin . In 
particular, on a plane surface normal to the propagation 
direction, 



and 



(:I n (p + S,t + T) :) = (: 7"(p,t) :), 



Pi 



(6) 
(7) 

P2 



where p is the transverse position vector, S 
and n = 1, 2, . . . 

For homogeneous and stationary fields described by 
positive nonsingular P distributions, the Schwarz in- 
equality implies that 



(T;I(p,t)I(p + 6,t + T) :> < (: I 2 (p,t) :), (8) 



that is, 



GW(S,t) <G^ 2 )(0,0). 



(9) 



Analogously to what was concluded from inequality (^), 
for fields that admit classical stochastic models, inequal- 
ity (j^) implies that photons are detected either spatially 
bunched or randomly spaced in a transverse detection 
screen. Violation of (g) implicates the possibility of quan- 
tum fields exhibiting spatial antibunching. Spatial anti- 
bunching of photons has been predicted by some authors 

@^!5iii2§. 



III. SPATIAL ANTIBUNCHING WITH 
DOWN-CONVERSION 

In this section we show that a field that violates in- 
equality (|^) can be generated by means of spontaneous 
parametric down-conversion. The experimental setup we 
are considering is shown in Fig. |[ A nonlinear birefrin- 
gent crystal is used to generate collinear entangled pho- 
ton pairs. The down-converted photons are the n incident 
on a bircfringent double-slit (see section III A ) and coin- 
cidences are detected by detectors D\ and D 2 . The pump 
beam is focused on the center of the plane of the double- 
slit, between the two slits. Interference filters are used 
such that the monochromatic approximation is valid. 

The following discussion refers to the basic geometry 
illustrated in Fig. [l], where a thin crystal is separated 
from an aperture plane by a distance s and the aperture 
plane is separated from a detection plane by a distance 
z. 

Using a treatment based in reference , in the parax- 
ial and monochromatic approximations, collinear SPDC 
generates a quantum state of the form |2q]: 




Aperture plane 



Detection plane 



FIG. 1: Illustration of the geometry, s is the crystal-aperture 
distance and z is the aperture-detector distance. 



with 



^Jdq 1 dq 2 $(q 1 ,q 2 )|q 1 ,cTi)|q 2 ,cT 2 ). (11) 



The coefficients C\ and C 2 are such that |C 2 1 \C\\. 
C 2 depends on the crystal length, the nonlinearity co- 
efficient, the magnitude of the pump field, among other 
factors. The kets |q-, cr,) represent Fock states labeled 
by the transverse component q ■ of the wave vector kj 
and the polarization aj of the down-converted photon 
j = 1,2. In this paper we consider type- II phase match- 
ing, in which case <j\ = e and a 2 = o where e (o) stands 
for extraordinary (ordinary) polarization. I^) is the two- 
photon component of the total quantum state. The func- 
tion ^(q!, q 2 ), which can be regarded as the normalized 
angular spectrum of the two-photon field |2q] , is given by 



« , (qi J q 2 ) 



2L 

TT 2 K 



u (Qi + cb) sinc 



^|qi - q 2 l 

4K 



'SPDC 



Ci|vac)+C 2 |V') 



(10) 



where v(q) is the normalized angular spectrum of the 
pump beam, L is the length of the nonlinear crystal in 
the z-direction, and K is the magnitude of the pump 
field wave vector. The integration domain D is, in prin- 
ciple, defined by the conditions q\ < kf and q 2 < 
However, in most experimental conditions, the domain 
in which <i>(q 1 ,q 2 ) is appreciable is much smaller than 
that. The state written above is not to be considered as 
a general expression for the SPDC process. Its validity is 
determined by experimental conditions, especially by the 
detection apparatus. As long as the monochromatic and 
paraxial approximations are valid, the results predicted 
by expression ( |ll| ) are in excellent agreement with ex- 
perience. Monochromatic approximation is guaranteed 
by the presence of narrow-band interference filters in the 
detection apertures, whereas paraxial approximation is 
guaranteed by keeping transverse detection regions much 
smaller than their distance from the crystal. 

We consider for now that the down-converted fields 
are incident on some sort of aperture, so as to produce 
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FIG. 2: Schematic diagram of spatial antibunching setup. 
An Ar laser pumps a BBO crystal, generating correlated pho- 
tons. The down-converted photons are incident on the bire- 
fringent double slit S and then the beamsplitter BS. The 
pump beam is focused on the double slit. Single and coinci- 
dence counts are registered with detectors Di and D2 . 



fourth-order interference in the absence of second-order 
interference. The reason for such a requirement is the 
following: Spatial photon antibunching is a fourth-order 
effect in a homogeneous field, that is to say, in a field that, 
according to |?j, does not show intensity patterns. With 
this scheme, we are seeking for a fourth-order interference 
pattern that depends only on x\ — X2, the relative position 
of detectors. Furthermore, this fourth-order interference 
pattern must have a minimum when x\ = x% in order 
to produce antibunching. Fourth-order spatial interfer- 
ence in the absence of second-order can be achieved in 
spontaneous parametric down-conversion by means of a 
double-slit whose slit separation is much greater than the 
transverse coherence length of the down-converted field, 
as reported by Fonseca et al. [Q. However, in refer- 
ence 1 26 the fourth-order correlation function, which is 
proportional to the coincidence rate, depends on x\ + X2 
instead of x\ — x%. In order to achieve a minimum of 
coincidences when x\ = X2, we have to introduce a phase 
difference of it between the two possibilities 
{photon 1 through slit 1, photon 2 through slit 2} and 
{photon 1 through slit 2, photon 2 through slit 1}. In 
our experiment, that phase difference was introduced by 
means of birefringent elements placed in front of each slit, 
as described later. After the aperture, the two-photon 
state can be written as 



MJ2 JJJJ d <li ^2 dqi dq' 2 *A (qi . qa) 
x^i.; (qi ~ Qi) ^4 (qj, - q 2 ) 



x|qi,0i)|q2,4) 



(13) 



where M is a normalization constant, $ / i(q 1 ,q 2 ) is the 
angular spectrum of the biphoton field on the aperture 



plane, that is, 



$ -A(qi,q 2 ) 



2L I L 

v(q t + q 2 ) sine ( — |q x - q 2 | 2 



TT 2 K 

x exp 



""■ ,+ *-srsl 



(14) 

T aa i (q) is the transfer function of the aperture, linking 
the incident field with transverse wave vector q and po- 
larization a with the scattered field with transverse wave 
vector q' and polarization a'. T aa i(q) is given by the 
Fourier transform of the aperture function A aa i(£). 

Since we are working with collincar SPDC with k% = 
fc 2 = \K, $a is written as 



<Mqi,q 2 ) 



2L 
n 2 K 



x exp 



u (qi + q 2 ) sm c (j^hi - q 2 
(|qi + q 2 | 2 + |qi-q 2 | 2 ) 



2K 



(15) 

where the irrelevant phase factor e %Ks is omitted. 

Using the orthonormal properties of the Fock states, 
we can define 



*(pi,p 2 ) = (v&c\V( P 2)<E)V(piM) 



(16) 



as the two photon coincidence detection amplitude, 
where 



V{p) 



„iki 



E 



dq a a (q) e a e l ^ p ~^ z) (17) 



is the monochromatic form of (Q) in the paraxial approx- 
imation and z is the distance between the aperture plane 
and the detection plane, as shown in Fig. [I. It is assumed 
that the polarization vector e is independent of q. The 
two-photon coincidence-detection probability for station- 
ary fields is proportional to the fourth-order correlation 
function with r = 0: 

V(pi,p 2 ) « G (2 ^( Pl , P2 ,0) - ||*( Pl ,p 2 )|| 2 . (18) 



A. The birefringent double-slit 

The birefringent double slit consists of two quarter- 
wave plates mounted in front of a typical double slit, 
such that each wave plate covers only one slit and their 
fast axes are orthogonal to one another, as shown in Fig. 
||. The slits are separated a distance d. With the plate- 
slit aperture oriented such that the slits and one fast axis 
are parallel to the e (y) direction and the other fast axis 
parallel to the o (x) direction, we can approximate the 
field-aperture functions by 

Aoo{£) = -i6(£ x -d/2) + 5(£ x + d/2) 
A ee {i) = 5(£ x - d/2) - i6(& + d/2) 

Aeo(Z) = 

A oe (S) = 0, (19) 



4 




Its angular spectrum is 



FIG. 3: The birefringent double slit. The quarter wave plates 
Qi and Q2 are aligned with orthogonal fast axes. S is a double 
slit with slits separation d. 



where £ x is the ir-component of £. The plate-slit aper- 
tures provide a controlled phase factor, that is, no phase 
will be added to a field with polarization parallel to the 
direction of the fast-axis of the wave-plate, while a field 
with perpendicular polarization will be modified by a 
phase factor of exp(— in/ 2). Thus, the phase factor de- 
pends on the polarization of the field as well as through 
which slit the field "passes" . 



B. The coincidence-detection probability amplitude 

Combining equations (|l3|-[r7|) and (|l^), we arrive at the 
following expression for coincidence-detection amplitude 
in the Fraunhofer approximation: 



* = * eo [e e <g> e G ] + * oe [e a <g> e e ], 



(20) 



where 



oc 



f \d kd 
{ cos -{qix + q2x) - + X V 



(21) 



where the "+" holds for \I> eo and the "-" holds for * oe . 

We assume that the pump field is a gaussian beam 
whose waist is located on the aperture plane: 



ua($) oc e 



v A {q) = v(q) exp 



is\K- 



2K 



oce-""' 2 / 4 , (23) 



where wo is the radius of the beam waist. Using (15) and 
in (|Tj) , it is straightforward to show that 



_f_<M 2 



kd. 

27 



(xx + x 2 ) 



=F sin 



kd 
27 



x 2 ) 



(24) 

It is interesting to note that the length L of the nonlin- 
ear crystal enters in the coincidence-detection amplitude 
only as a multiplicative constant. It is clear from expres- 
sion ( p4[ ) above, that the fulfillment of the homogeneity 
condition urn for n — 2 in the x-direction depends on the 

factor e y2 ™a> . If wq d, the dependence on X\ + X% 
disappears and transverse field on the detection plane 
can be considered as homogeneous. This is the reason 
why the pump beam must be focused on the center of 
the double slit. In this case, 



*eo(Pl,P2) = -*oe(Pl,P2) OC Sm 



kd 

Yz {xi 



x 2 ) 



Thus, the coincidence detection probability is 

'kd. 



V(px,P2) oc 1 - cos 



-{xi - x 2 ) 



(25) 



(26) 



(22) 



When x\ — x 2 , the coincidence count rate is zero 
and increases with X\ — x 2 until [x,\ — x 2 )kd/z = 
±7r/2. Therefore, the fourth-order correlation function 
G^ 2 ' 2 '(pi, P%i t), which is proportional to the coincidence 
detection probability V(pi, p 2 ) 7 does not have a maxi- 
mum at xi — x 2 . This contradicts thus characteriz- 
ing spatial antibunching of photons. 



IV. DISCUSSION AND CONCLUSION 

We have shown the theoretical background behind the 
spatial antibunching of photons using parametric down- 
conversion. It may be instructive for the reader to 
compare the experiment analyzed here with its classi- 
cal counterpart. In this context, the single count de- 
tection rate R c i{x) should be proportional to the clas- 
sical average intensity (I(x)), whereas the coincidence 
count rate C c i(xi,x 2 ) should be proportional to the 
intensity-intensity (or the fourth-order) correlation func- 
tion (I(xi)I(x 2 )) . The single count detection rate of 
down-converted light in the presence of a double-slit has 
been studied in previous works |p6| , p7| , p8f . In refer- 
ence p7| it was demonstrated that in terms of its single 
count rate, SPDC behaves like a classical Schcll-modcl 
extended light source. In our experiment, the transverse 
coherence length being shorter than the slits separation 
and shorter than the slits widths themselves, the single 
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count rate is given by the classical expression for incoher- 
ent illumination, which can be approximated by a gaus- 



R c i{x) cx e 



(27) 



where a 



ka 



and a is the width of the slits. Since 



the transverse detection range x max — x m i n is much 
shorter than the width of this gaussian profile for the 
slits-detectors distance considered, the single count rate 
is fairly constant over de detection range Q. By an- 
other side, the coincidence detection rate due to a clas- 
sical source is totally different from the observed with 
down-converted light. Perhaps, the best classical model 
for type-II SPDC is a superposition of two extended light 
sources orthogonally polarized and correlated in inten- 
sity. After the light is diffracted by the birefringent 
double-slit, the calculation of the classical fourth-order 
correlation function is quite similar to the case of the 
Hanbury Brown - Twiss intensity interferometer p9| . 
Classical intensity intcrfcrometry is known to be insensi- 
tive to phase. Therefore, the birefringent elements have 
no effect on the predicted fourth-order correlation, that 



C c i (xi , x%) oc 1 + v cos 



kd 

— (Xl - x 2 ) 



(28) 



The visibility v is in the range < v < ^ and depends 
on the statistics of the source. It is clear from expression 
( p8| ) above that C c i(xi, x 2 ) predicts spatial bunching, as 
expected from any classical light source. In view of the 
above analysis, the results presented here describe an en- 
tirely quantum fourth-order interference effect, with no 
classical analog |3(J . In addition to rendering further in- 
terest in the study of non-classical states of light, spatial 
antibunching promises to be a useful tool in quantum 
imaging and quantum information technologies. 
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